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Abstract. - We present a computer simulation model of Wheeler's delayed choice experiment 
that is a one-to-one copy of an experiment reported recently (V. Jacques et ai., Science 315, 966 
(2007)). The model is solely based on experimental facts, satisfies Einstein's criterion of local 
causality and does not rely on any concept of quantum theory. Nevertheless, the simulation 
model reproduces the averages as obtained from the quantum theoretical description of Wheeler's 
delayed choice experiment. Our results prove that it is possible to give a particle-only description 
of Wheeler's delayed choice experiment which reproduces the averages calculated from quantum 
theory and which does not defy common sense. 



Introduction. — According to the wave-particle du- 
ality, a concept of quantum theory, photons exhibit both 
wave and particle behavior depending upon the circum- 
stances of the experiment [1]. In 1978, Wheeler proposed 
a gedanken experiment [2], a variation on Young's double 
slit experiment, in which the decision to observe wave or 
particle behavior is made after the photon has passed the 
slits. The pictorial description of this experiment defies 
common sense: The behavior of the photon in the past is 
said to be influenced changing the representation of the 
photon from a particle to a wave. 

Recently, Jacques et ai. reported an almost ideal ex- 
perimental realization of Wheeler's delayed choice exper- 
iment [3]. The experimental set-up (see Fig. 1) consists 
of a single-photon source, a Mach-Zehnder interferome- 
ter, with at the output side a beam splitter (BS outpu t) of 
which the presence can be controlled by a voltage applied 
to an electro-optic modulator (EOM) and detectors [3]. 
Although the detection events are the only experimental 
facts, the pictorial description of Jacques et ai [3] is as fol- 
lows: The decision to apply a voltage to the EOM is made 
after the photon has passed BSi„ ptl t but before the pho- 
ton enters BS output- If no voltage is applied to the EOM 
(open configuration), then the arrival of a photon at ei- 
ther detector clearly gives which-way information about 
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Fig. 1: Schematic diagram of the experimental setup for 
Wheeler's delayed-choice gedanken experiment [3]. PBS: Polar- 
izing beam splitter; HWP: Half-wave plate; EOM: electro-optic 
modulator; RNG: Random number generator; WP: Wollaston 
prism; P,S: Polarization state of the photon; Do, D%: Detec- 
tors. 



the photon within the interferometer (particle behavior), 
with 50% arriving from either path. When the voltage is 
applied (closed configuration), the paths interfere and it 
is impossible to know which path the photon took (wave 
behavior). Accordingly, the detectors register an interfer- 
ence pattern. 

The outcome of delayed-choice experiments [3 8], that 
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is the average results of many detection events, is in agree- 
ment with quantum theory. However, the pictorial de- 
scription [3] defies common sense: The decision to apply 
a voltage to the EOM after the photon left BS input but 
before it passes BS output , influences the behavior of the 
photon in the past and changes the representation of the 
photon from a particle to a wave [3]. On the other hand, 
quantum theory does not describe single events [1]. There- 
fore, it should not be a surprise that the application of 
concepts of quantum theory to the domain of individual 
events may lead to conclusions that are at odds with com- 
mon sense. 

In this Letter, we describe a model that, when imple- 
mented as a computer program, performs an event-by- 
event simulation of Wheeler's delaycd-choice experiment. 
Every essential component of the laboratory experiment 
(PBS, EOM, HWP, Wollaston prism, detector) has a coun- 
terpart in the algorithm. The data is analyzed by counting 
detection events, just like in the experiment [3]. The simu- 
lation model is solely based on experimental facts, satisfies 
Einstein's criterion of local causality and does not rely on 
any concept of quantum theory or of probability theory. 
Nevertheless, our simulation model reproduces the aver- 
ages obtained from the quantum theoretical description of 
Wheeler's delayed choice experiment but as our approach 
does not rely on concepts of quantum theory and gives a 
description on the level of individual events, it provides a 
description of the experimental facts that does not defy 
common sense. In a pictorial description of our simula- 
tion model, we may speak about "photons" generating 
the detection events. However, these so-called photons, as 
we will call them in the sequel, are elements of a model 
or theory for the real laboratory experiment only. The 
experimental facts are the settings of the various appara- 
tuses and the detection events. What happens in between 
activating the source and the registration of the detec- 
tion events is not measured and is therefore not known. 
Although the photons "know" exactly which route they 
followed in the closed configuration of the interferometer 
(we can always track them during the simulation), they 
build up an interference pattern at the detector. The 
appearance of an interference pattern is commonly con- 
sidered to be characteristic for a wave. In this Letter, we 
demonstrate that, as in experiment, it can also be build up 
by many photons. These photons have which- way infor- 
mation, never directly communicate with each other and 
arrive one by one at a detector. 

To head off possible misunderstandings, the work pre- 
sented here is not concerned with the interpretation or an 
extension of quantum theory. We adopt the point of view 
that quantum theory has nothing to say about individual 
events [1]. The fact that there exist simulation algorithms 
that reproduce the results of quantum theory has no direct 
implications to the foundations of quantum theory: These 
algorithms describe the process of generating events at 
a level of detail that is outside the scope of what current 
quantum theory can describe. The event-based simulation 



approach that we describe in this Letter is unconventional 
in that it does not require knowledge of the probability dis- 
tribution obtained by solving the quantum problem. The 
averages given by quantum theory are obtained through a 
simulation of locally causal, classical dynamical systems. 
The key point of these dynamical systems is that they are 
built from units that are adaptive. 

It is common practice to use the framework of Kol- 
mogorov 's probability theory to construct probabilistic 
models of phenomena that cannot (yet) be described by 
a deductive theory. Although Kolmogorov 's probability 
theory provides a rigorous framework to formulate such 
models, there are ample examples that illustrate how easy 
it is to make plausible assumptions that create all kinds of 
paradoxes, also for every-day problems that have no bear- 
ing on quantum theory at all [9-13]. Subtle mistakes such 
as dropping (some of the) conditions, mixing up the mean- 
ing of physical and statistical independence, and changing 
one probability space for another during the cause of an ar- 
gument can give rise to all kinds of paradoxes [14-16]. To 
avoid these potential pitfalls, in our simulation approach 
we strictly stay in the domain of integer arithmetic, that 
is we do not invoke any concept of probability theory. 

This Letter builds on earlier work [17-25] that demon- 
strated that it may be possible to simulate quantum phe- 
nomena on the level of individual events without invoking 
a single concept of quantum theory. Specifically, we have 
demonstrated that locally-connected networks of process- 
ing units with a primitive learning capability can simu- 
late event-by-event, the single-photon beam splitter and 
Mach-Zehnder interferometer experiments of Grangier et 
al. [26] and Einstein-Podolsky-Roscn experiments with 
photons [27-29]. Furthermore, we have shown that this 
approach can be generalized to simulate universal quan- 
tum computation by an event-by-event process [19,20]. 
The algorithms used in our earlier work [17-25] cannot be 
used to simulate Wheeler's delayed choice experiment [3]. 
The latter uses components that respond to the polariza- 
tion and/or the difference in path length and the algo- 
rithms used in our earlier work cannot handle both fea- 
tures simultaneously. In contrast, the more general algo- 
rithms described in this Letter can also simulate all the 
experiments covered in our earlier work. Our cvent-by- 
event simulation approach rigorously satisfies Einstein's 
criterion of local causality and builds up the final out- 
come that agrees with quantum theory event-by-event, as 
observed in real experiments. 

Simulation model. — The simulation algorithm can 
be viewed as a message-processing and message-passing 
process: It routes messengers through a network of units 
that process messages. In a pictorial description, the pho- 
ton is the messenger, carrying a message representing its 
phase and polarization. The processing units play the 
role of the components of the laboratory experiment and 
the network represents the complete experimental setup. 
Some processing units consist of an input stage (a stan- 
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dard linear adaptive filter), a transformation stage and 
an output stage. The input (output) stage may have sev- 
eral channels at (through) which messengers arrive (leave). 
Other processing units are simpler in the sense that the 
input stage is not necessary for the proper functioning of 
the device. A message is represented by a vector. As a 
messenger arrives at an input channel of a processing unit, 
the input stage updates its internal state, represented by 
a vector, and sends the message together with its internal 
state to the transformation stage that implements the op- 
eration of the particular device. Then, a new message is 
sent to the output stage, using a pseudo-random number 
to select the output channel through which the messenger 
will leave the unit. We use pseudo-random numbers to 
mimic the apparent unpredictability of the experimental 
data only. The use of pseudo-random numbers is merely 
convenient, not essential. At any given time, there is only 
one messenger being routed through the whole network. 
There is no direct communication between the messengers. 

In the experimental realization of Wheeler's delayed 
choice experiment by Jacques et al. [3] linearly polarized 
single photons are sent through a polarizing beam splitter 
(PBS) that together with a second, movable PBS forms 
an interferometer (see Fig. 1). The network of processing 
units is a one-to-one image of the experimental setup [3] 
and is therefore not shown. We now describe each of the 
components of the network in detail. 

Messenger. Photons are regarded as messen- 
gers. Each messenger carries a message repre- 
sented by a six-dimensional unit vector yk, n = 
( cos V^„, sinV>£„, cos V^„, sin^„, cos£ fe ,„, sin£ fei „). The 
superscript H (V) refers to the horizontal (vertical) compo- 
nent of the polarization and ipj? n , n-> ano - represent 
the phases and polarization of the photon, respectively. 
It is evident that the representation used here maps one- 
to-one to the plane-wave description of a classical electro- 
magnetic field [30] , except that we assign these properties 
to each individual message, not to a wave. The subscript 
n > numbers the consecutive messages and k = 0, 1 la- 
bels the channel of the PBS at which the message arrives 
(see below). 

Polarizing beam splitter. The processor that performs 
the event-by-event simulation of a PBS is depicted in 
Fig. 2. It consists of an input stage, a simple deterministic 
learning machine (DLM) [17-20], a transformation stage 
(T), an output stage (O) and has two input and two out- 
put channels labeled with k = 0,1. We now define the 
operation of each stage explicitly. 
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Fig. 2: Diagram of a DLM-based processing unit that performs 
an event-based simulation of a polarizing beam splitter (PBS). 
The solid lines represent the input and output channels of the 
PBS. The presence of a message is indicated by an arrow on 
the corresponding channel line. The dashed lines indicate the 
data flow within the PBS. 
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ternal vector x„ = (xo in , x\ >n ), where xo, n + x\^ n = 1 
and Xk, n > for k = 0, 1 and all n. These seven two- 
dimensional vectors are labeled by the message num- 
ber n because their contents are updated every time 
the DLM receives a message. Note that the DLM 
stores information about the last message only. The 
information carried by earlier messages is overwritten 
by updating the internal registers. 

Upon receiving the (n + l)th input event, the DLM 
performs the following steps: It stores the first two 
elements of message yk,n+i in its internal register 
Y "n+i = ( C k,n+i> S k,n+i)> thc middle two elements 
of y fe ,„+i in Y^„ +1 = (Cl n+1 ,Sl n+1 ), and the last 
two elements of y fc ,„ + i in Y£„ +1 = (C£ n+1 , S£ n+1 ). 
Then, it updates its internal vector according to the 
rule [17] 



Ci, n +i = axi,„ + (1 - a)d it k, 



(1) 



where < a < 1 is a parameter that controls the 
learning process [17]. Note that by construction 
Zo,n+i + xx.n+i = 1, x , n+ i > and x^ n+ i > 0. 



Input stage: The DLM receives a message on either 
input channel or 1, never on both channels simulta- 
neously. The arrival of a message on channel (1) is 
named a (1) event. The input events are represented 
by the vectors v„ = (1,0) or v„ = (0, 1) if the nth 
event occurred on channel or 1, respectively. The 
DLM has six internal registers Y^ n = (Cj? n , Sf? n ), 



Transformation stage: Thc second stage (T) of 
the DLM-based processor takes as input the data 
stored in the six internal registers Y^ n+1 = 
(n H q h \ v v — in v q v \ 

\^k,n+l' °k,n+l>' 1 fe,n+l — ^/c^+l ' °k,n+l )' 

Y k,n+i = ( C k,n+1> S k,n+i) and in thc internal vector 
= {xo. ,n+i) 2-1,71+1 ) 9-nd combines the data into 
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the eight-dimensional vector 



T = 



, 5'o.n+lCo,ri+l\/ a;:0 '™+ 1 

_ qV cP i 

D l,ri+1 J l,n+1 \/ x l,n+l 

^l^n+l ^l^n+l y/ x l.n+l 

~^X,n+l^Q,n+l\/ x 0.n+l 
V ^0,n+l^0,n+W X 0-n+l / 



(2) 



Rewriting the vector T as 



T = 



^ ipO,n+l + ^ofn+l) < ^a^i+l\/ x 0.n+l ^ 

+ ^lfn+l) ^l,n+l y'Xi.n+i 
V * (Qvi+1 + ^n+l) ^ari+l\/^0^+l / 



of 



(3) 



shows that the operation performed by the transfor- 
mation stage T corresponds to the matrix- vector mul- 
tiplication in the quantum theoretical description of 
a PBS, namely 
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where (aff , o y , af , a\) and (bff , 6 y , 6f , b\) denote 
the input and output amplitudes of the photons with 
polarization H and V in the and 1 channels of a 
PBS, respectively. Note that in the quantum optical 
description of a (polarizing) beam splitter the vac- 
uum field must be included. In our simulation model, 
there is no need to introduce the concept of a vacuum 
field. 

Output stage: The final stage (O) sends the message 



w 



where 

W ,n+1 
Wl,n+1 
W2,n+1 
W3,n+1 

S 0,n+1 



( W , n +l/S0,n+l \ 
Wl,n+l/ S0,n+1 
W2,n+l/si, n +l 
^3,7l+l/ Sl.n+1 
,n+l 

\ Sl,n+l/S2,n+l / 



Co,n+lCo,n+l\/ X 0,n+l, 
$0,71+1 Co.n+1 \/ x 0,n+l, 



(5) 
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2,rt+l "T" w 3,n+l 



s l,n+l — 

g 2,n+l = \J wl n+1 + W 2 ln+1 + wl n+1 



J 3,n+1' 

(6) 



(7) 



through output channel if Wq, 71+1+^1,71+ 1 > ^ where 
< r < 1 is a uniform pseudo-random number. Oth- 
erwise, if Wq „ +1 < r, the output stage sends 
through output channel 1 the message 

/ ZO,n+l/to,n+l \ 
Zl.n+l/to.n+1 
Z2,n+l/h,n+l 
Z3,n+l/tl,n+l 

^0,71+1/^2,71+1 

\ *l,n+l/*2,n+l / 

where 

^o.ti+i = C^ n+1 C[ n+lv /xi^ n+ i, 

Z2,n+1 — ~So,n+lSo,n+l\/ X 0,n+l, 
^3,71+1 = C'0,n+1^0,n+l\ /x 0,n+l, 



^0,71+1 
^1,71+1 
^2,71+1 



),7l+l ^ z l,n+lt 



'2,71+1 



+ 2 



3,n+lJ 



'0,71+1 



-l,n+l 



2,71+1 



+ 2 



3,71+1- 



(8) 



As mentioned earlier, the use of pseudo-random numbers 
to select the output channel is not essential. We use 
pseudo-random numbers to mimic the apparent unpre- 
dictability of the experimental data only. Instead of a 
uniform pseudo-random number generator, any algorithm 
that selects the output channel in a systematic manner 
might be employed as well. This will change the order in 
which messages are being processed but the content of the 
messages will be left intact and the resulting averages do 
not change significantly. 

Half wave plate (HWP). This device performs a rota- 
tion of the polarization of the photon [30]. If the polariza- 
tion of an incoming photon is at an angle 9 with respect to 
the optical axis of the HWP then, after passing the HWP, 
the polarization of the photon has been rotated by an an- 
gle 29. In order to change horizontal polarization into 
vertical polarization, or vice versa, a HWP is used with 
its optical axis oriented at ir/4 . The HWP does not only 
change the polarization of the photon, but also its phase 
as can be seen from the wave mechanical description 



b H 
b v 
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71 



cos 26» sin 29 
sin 29 -cos 26* 



(9) 



1,71+1' 



As a result, for the case 6 = 7r/4 , the polarization of the 
photon is rotated by an angle tt/2 and its phase is changed 
by -tt/2. 
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Fig. 3: Simulation results generated by the DLM net- 
work of the PBS shown in Fig. 2. Input chan- 
nel receives (cos ipQ , sin ipo , cos 4>o , sin^J', cos 8, sin 9) — 
(1, 0, 1, 0, cos 8, sin 8). Input channel 1 receives no events. After 
each set of N = 10000 events, 8 is increased by 15°. Squares 
and circles give the simulation results for the normalized in- 
tensities No/N and N\/N as a function of 9. Dashed lines 
represent the results of quantum theory. 



Electro-optic modulator (EOM). This device rotates 
the polarization of the photon by an angle depending on 
the voltage applied to the modulator. In the laboratory 
experiment [3], the EOM is operated such that when a 
voltage is applied the EOM acts as a HWP that rotates the 
input polarizations by 7r/4. In the simulation a pseudo- 
random number is used to decide to apply a voltage to the 
EOM or not. Also here we use a pseudo-random number to 
mimic the experimental procedure to control the EOM [3] . 
Any other (systematic) sequence to control the EOM can 
be used as well. 

Wollaston prism. This device is a PBS with one input 
channel and two output channels and is simulated as the 
PBS described earlier. Messengers arrive at one and the 
same input channel only. 

Detection and data analysis procedure. Detector Dq 
(Di) registers the output events at channel (1). During 
a run of N events, the algorithm generates the data set 

T = {x n , A n \n = 1, .... JV; (j> = fa - M , (10) 

where x n = 0,1 indicates which detector fired (Dq or Di), 
and A n = 0, 1 is a pseudo-random number that is chosen 
after the nth message (=photon) has passed the first PBS, 
determining whether or not a voltage is applied to the 
EOM (hence whether the configuration is open or closed) . 
The angle 4> denotes the phase shift between the two inter- 
ferometer arms. This phase shift is varied by applying a 
plane rotation on the phase of the particles entering chan- 
nel of the second PBS. This corresponds to tilting the 
second PBS in the laboratory experiment [3]. For each 
phase shift <fi and interferometer configuration (open or 



Fig. 4: Event-by-event simulation results of the experimental 
realization of Wheeler's delayed-choice gedanken experiment. 
Open (closed) markers correspond to data for the open (closed) 
configuration of the interferometer. Squares and circles give 
the simulation results for the normalized intensities iVo /N and 
N\/N as a function of the phase shift cj> , No (Ni) denoting 
the number of events registered at detector Do (Di). For each 
value of <f>, the number of input events N = 10000. The number 
of detection events per data point is approximately the same as 
in experiment. The simulation data is in qualitative agreement 
with experiment: See Fig. 3 of Ref . [3]. Dashed lines represent 
the results of quantum theory (Malus law). 



closed) the number of (1) output events No (Ni) 
culated. 



is cal- 



Simulation results. — The algorithm described 
above directly translates into a simple computer program. 
For a fixed set of input parameters, each simulation takes 
a few seconds on a present-day PC. All simulations are 
carried out with a = 0.99 [17]. 

We first demonstrate that our model for the PBS re- 
produces Malus' law. In this simulation, we send mes- 
sengers to one input channel, say channel 0, only. This 
implies that the registers that are connected to chan- 
nel 1 will not change during the simulation. Figure 3 
shows a representative set of event-based simulation re- 
sults for the PBS modeled by the processor depicted in 
Fig. 2. The data set is produced with input messages 
(cos ipQ , sin tpQ , cos tp^ > sm > cos £o , sin £ ) • The values 
of tpQ and ifo are fixed but irrelevant otherwise. Also the 
value of £o is irrelevant. It is clear that the intensities in 
both output channels obey Malus' law. 

Next, we build a network that contains all the optical 
components of the laboratory experiment [3] (see Fig. 1). 
Before the simulation starts we set Xo = (2:0,0, ^1,0) = 
(r, 1 — r), where r is a uniform pseudo-random number. In 
a similar way we use pseudo-random numbers to initialize 

Y o% Y ^o> Y £o> Y fo, Y i.o and Y i,o- In tllis simulation, 
we send messengers to one input channel (see Fig. 1). The 
HWP in BS ou tput changes the phases and also interchanges 
the roles of channels and 1. Disregarding a few excep- 
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tional events, the PBS in BS OM ij,„t generates messages in 
one of the channels only. 

Representative results of an event-by-event simulation 
of this network, a one-to-one image of Wheeler's delayed 
choice experiment [3], are shown in Fig. 4. The simulation 
data are in quantitative agreement with the averages cal- 
culated from quantum theory and in qualitative agreement 
with experiment [3]. 

Conclusion. — In this Letter, we have proven that it 
is possible to give a particle-only description for both the 
open and closed interferometer configuration in Wheeler's 
delayed choice experiment that (1) reproduces the aver- 
ages calculated from quantum theory, (2) satisfies Ein- 
stein's criteria of realism and local causality, (3) does not 
rely on any concept of quantum theory or of probability 
theory, and (4) is not in conflict with common sense. 
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